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1 Introduction 

Elliptic functions are known to appear in many problems, applied and the- 
oretical. However, a lesser known application is in the study of exact solu- 
tions to Einstein's gravitational field equations in a Friedmann-Robertson- 
' Cj': Lemaitre- Walker (FRLW) cosmology [H H El El [151 [16] . We will show ex- 

^1 plicitly how Jacobi and Weierstrass elliptic functions arise in this context, 

^1 and will additionally show connections with theta functions. In section 2, 

we review the definitions of various elliptic functions. In section 3, we record 
relations between elliptic functions and theta functions. In section 4 we in- 
troduce the FRLW cosmological model and we then proceed to show how 
QQ I elliptic functions appear as solutions to Einstein's gravitational equations in 

^ I sections 5 and 6. The author thanks Floyd Williams for helpful discussions. 

00 ^ 

^ 2 Elliptic functions 

o ■ 

An elliptic integral is one of the form J R{x, ^J P{x))dx where P{x) is a 
I polynomial in x of degree three or four and i? is a rational function of its 

I arguments. Such integrals are called elliptic since an integral of this kind 

arises in the computation of the arclength of an ellipse. Legendre showed 
that any elliptic integral can be written in terms of the three fundamental or 
normal elliptic integrals 



F{x, k) 



dt 



V(l-t2)(l-A;2t2) 



(2.1) E{x,k)''d- I \r-Y:^dt 



1 - A;2t2 



1 



Il{x, o? ^ k) 



dt 



which are referred to as normal eUiptic integrals of the first, second and third 
kind, respectively. The parameter a is any real number and k is referred to 
as the modulus. For many problems in which real quantities are desired, 
< k,x < 1, although this is not required in the above definitions (for A; = 
and k — 1, the integral can be expressed in terms of elementary functions 
and is therefore pseudo-elliptic) . 

EUiptic functions are inverse functions of elliptic integrals. They are 
known to be the simplest of non-elementary functions and have applications 
in the study of classical equations of motion of various systems in physics 
including the pendulum. One can easily show that if f{u) denotes the in- 

verse function of an elliptic integral ^ ^))d., then since 

y{f{u)) = u^y'{f{u)) = l/r{u), 

(2.2) fiuf = 



R{f{u),,/P{fWW' 

The Jacobi elliptic function sn{u, k) is the inverse of F{x, k) defined 
above, and eleven other Jacobi elliptic functions can be written in terms 
of sn{u, k): cn{u, k) and (in(ti, k) satisfy 

(2.3) sn^u -\- cn^u — 1 and k'^sn^u + dn^u — 1 

respectively, and 



(2.4) 



By (2.1), (2.2), (2.3) and (2.4) one can see that each one of the Jacobi elliptic 
functions satisfy the differential equation 

(2.5) f'^uf + af{uf + hf{uY = c 

for some a, 6, c in terms of the modulus k according to the following table. 
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1 - 2P 
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nc{u, k) 


1-2P 




-k' 


nd{u, k) 
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-1 


sc{u, k) 
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sd{u, k) 


1-2^ 


e{i-k^) 
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cd{u, k) 


i + e 
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cs{u, k) 


e-2 


-1 


1 - p 


ds{u, k) 


1-2^ 


-1 


A;2(A;2 - 1) 


dc{u, k) 


1 + e 
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Note that the a- value for f{u) = dn{u, k) seen above corrects an error in [6]. 
The Weierstrass elliptic function 



(2.6) p{z; uJi,uj2) 



E 



{m,n)eZxZ-{(0,0)} 



{z — muji — nuj2Y {rnuji + nuj2y 



is a doubly-periodic elliptic function of ^ e C with periods uji,uj2 G C such 
that lm{uj\lijj2) > 0. p is the inverse of the elliptic integral 



(2.7) 



P \x;g2,g3) 



v/4t3 - g^t - g-i 



dt 



where 5'2,5'3 € C are known as Weierstrass invariants. Given periods a;i,a;2, 
the invariants are 



92 



93 



60 Y, 



1 



(m,n)eZxZ-{(0,0)} 



140 



{muJi + nuj2y- 
1 



(m,n)eZxZ-{{0,0)} 

Alternately given invariants 92, gs, periods uji,uj2 can be constructed if the 
discriminant A g2—21g1 is nonzero - that is, when the Weierstrass cubic 
4:t^—92t—93 does not have repeated roots (see 21-73 of Whittaker and Watson 
). In this lecture we refer to p{z; Ui, U2) as either p{z) or p{z; g2, g^) and 
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consider only cases where the invariants are such that 7^ '^'^dl- By (2.2) 
and (2.7) one can see that the Weierstrass elhptic function satisfies 

(2.9) p'{zf = Ap{zf-g^p{z)~g:,. 

Note that in Michael Tuite's lecture in this volume, u}m,n there is equal to 
muJi + nu2 here with our uji,uj2 specialized to 27rzr and 27rz in his lecture. 

In the special case that the discriminant A > 0, the roots of the Weier- 
strass cubic 4:t^—g2t — g3 are real and distinct, and are conventionally notated 
by ei > 62 > 63 for ei + 62 + 63 = (since the cubic contains no squared 
term). In this case At^ — g2t — g^ = 4(t — ei){t — e2){t — 63) so that the 
Weierstrass invariants are given in terms of the roots by 

(2.10) 6-2 = -4(eie2 + 6163 + 6263), 5'3 = 4616263. 

For 61 < z G R, p can then be written in terms of the Jacobi elliptic functions 
in any one of the following ways 

piz;g2,g3) = 63 + 7^ns^(72;, /c) 

(2.11) p{z;g2,g3) = e2 + l'^ds'^{jz,k) 

Piz;g2,g3) = 61 +7^cs^(7Z, fc) 

where 7^ 61 — 63 and the modulus k is such that k"^ = f^5f| (similar 
equations hold if z G R is in a different range in relation to the real roots 
61,62,63, and alternate relations hold for non-real roots when A < 0, see 
chapter II of Greenhill |12]). 

Note that the Jacobi elliptic functions solve a differential equation which 
contains only even powers of f{u), and p solves an equation with no squared 
or quartic powers of p. Weierstrass elliptic functions have the advantage of 
being easily implemented in the case that the cubic 4:t^—g2t—g3 is not factored 
in terms of its roots. Elliptic integrals of type J R{x) / a/ P{x)dx, where P(x) 
is a cubic polynomial and i? is a rational function of x, can be written in 
terms of three fundamental Weierstrassian normal elliptic integrals although 
we will not record the details here (see Appendix of Byrd and Friedman [8]). 
In section 6 we will see a method which allows one to write the elliptic integral 

1/ \jF{t)dt, for F{t) a quartic polynomial, in terms of the Weierstrassian 
normal elliptic integral of the first kind (2.7) by reducing the quartic to a 
cubic. 
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3 Jacobi theta functions 



Jacobi theta functions are functions of two arguments, z G C a complex 
number and r G H in the upper-half plane. Every elliptic function can be 
written as the ratio of two theta functions. Doing so elucidates the mero- 
morphic nature of elliptic functions and is useful in the numerical evaluation 
of elliptic functions. One must be cautious with the notation of theta func- 
tions, since many different conventions are used. We will use the notation of 
Whittaker and Watson 1201 to define 



(3.1) 



Oi{z,t) 



03{z,t) 



def. 



def. 



def. 



def. 



2q^/^ ^(-l)"g"("+^)szn((2n + l)z) 

n=0 

oo 

2g^/^ g"("+i)cos((2n + l)z) 

n=0 
oo 

1 + 2 J]g"'cos(2n;2) 

n=l 

oo 

1 + 2^(-l)V'cos(2n2). 



n=l 



def. 



where q =' e'^*^ is called the nome. We also define the special values 6', 
In terms of theta functions, the Jacobi elliptic functions are 



def. 



(3.2) 



sn{u, k) 
cn{u, k) 
dn{u, k) 



03 




O2 


e,{u/eiT) 


O4 


e,{u/elr) 


02 




Oa 


d3{u/elT) 


O3 





where r is chosen such that k"^ = ^'2/^3. By 22 ■ 11 of Whittaker and Watson 
[20] . if < fc^ < 1, there exists a value of r for which the quotient 6'|/6'| = k"^. 



5 



4 The FRLW cosmological model 



The Friedmann-Robertson-Lemaitre- Walker cosmological model assumes 
that our current expanding universe is on large scales homogeneous and 
isotropic. On a d + 1— dimensional spacetime this assumption translates 
into a metric of the form 

(4.1) ds' = -df + a{tf f + r^ciOLi 

where a{t) is the cosmic scale factor and k' e {—1,0,1} is the curvature 
parameter. 

Einstein field equations Gij = —KdTij + Agij then govern the evolution 

of the universe over time. In these equations, the Einstein tensor Gij 
Rij — \Rgij is computed directly from the metric gij by calculating the Ricci 
tensor Rij and the scalar curvature R. Also — 87rGd where Ga is a 
generalization of Newton's constant to d + 1— dimensional spacetime and 
A > is the cosmological constant. The form of the energy-momentum 
tensor Tij depends on what sort of matter content one is assuming, and in 
this lecture will be that of a perfect fluid - that is, Tij = (p + p)giogjo + pgij 
where p{t) and p{t) are the density and pressure of the fluid respectively. 

For the metric (4.1) with a perfect fluid energy-momentum tensor, Ein- 
stein's equations are 

def ^ 

for B.{t) = a{t)/a{t) and where dot denotes differentiation with respect to 
t. In this lecture, only equation (i) will be required to relate the cosmic scale 
factor a(t) to elliptic and theta functions. 

We begin by rewriting equation (i) in terms of conformal time rj by defin- 
ing the conformal scale factor a{r]) a{f{ri)) where f{ri) is the inverse 
function of ri{t) which satisfies f]{t) — l/a{t). In terms of a{ri), (i) becomes 

(4.2) a'ir^f = Aa{vY + ^dP{f{v)HvY - k'a{r^f 

~ def. ~ def . 

where we use notation A =' 2K/d{d — 1), — 2Kd/d{d — 1) and we take 
spatial dimension d> 1. 
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5 FRLW and Jacobi elliptic and theta 
functions 

In general, if f{u) is a solution to f'{uf + af{uy + hf{uY — c, then 
g{u) = Pf{au) is a solution to 

(5.1) g'^uf + Ag{uf + Bg{uY = ^ 

an 

for a — and (3 — where we may choose either the positive or 

negative square root for each of a and p. We will construct solutions to 
the conformal time FRLW Einstein equation (4.2), given that Jacobi elliptic 
functions solve (2.5), and also proceed to write these solutions in terms of 
theta functions by the relations in (3.2). 

For the special case of density p{t) = with D > 0, (4.2) becomes 

(5.2) a'{r])^ + k'a{r]f - la{r]Y = KdD. 

We construct real-valued solutions to this equation as above in terms of Ja- 
cobi elliptic functions, by restricting to special values of k' and D — J^^~^ > 
for which the ratios a/ A and b/B are positive. 

In the case of positive curvature k' — 1 and D — — zr^ for any 

;j5A(i+jfc2)2 

< A; < 1, the conformal time equation (5.2) becomes 

(5.3) a\rif + a{rif - Aa{ri)^ = ^^-^^ 

A (1 + A;2)2 

so that A = 1 and B = —A. By considering Jacobi elliptic functions for 
which ^ = ^ and -f = are positive, some solutions to (5.3) in terms of 
Jacobi functions (and their equivalent theta function representations) are 

(5.4) 

asniv) = I sn A; ' - 



A(i + A;2) Vv/T+F' y JjJ^l^) o^{r]/./efTel,T) 



ansiv) = I 'ns . k 



adciv) = I dc' ' - 



^A(l + k^) Vv^TTF' ; ^A(^34 + ei) Hv/Voi + oir) 

where r is chosen such that A;^ = ^|/^|. 

The first two solutions, agniv) and ans{Tj), reduce to hypcrbohc trigono- 
metric functions in the case of modulus k = 1, since sn{u, 1) = tanh{u) and 
ns{u, 1) = coth{u). That is, two additional solutions in terms of elementary 
functions are 

(5.5) ai{r]) = ^= tanh{r]/V2) 

V2A 

02(77) = ^= coth{rj/\/2). 
V2A 

For these two solutions, one may solve the differential equation f]{t) = 
l/a{t) — 1/0(77 (i)) for v{t): a-iid therefore obtain the cosmic scale factor 
a(t) = a{rj{t)) which solves the Einstein field equation (i) in section 4 for 
p{t) = D/a{tY with special values k' — 1 and D — k'^/AKd{l + k'^Y- Doing 
so. we obtain 



(5.6) 




a\{t) = ai{r){t)) = tank ( In 



'2A 

02 (i) = a2{r]{t)) — coth I In 



'2A 

for t > 0. 

In the case of negative curvature k' — —1 and D — for any 

< k < 1, the conformal time equation (5.2) becomes 

(5.7) a'irif - 0(77)^ - Aa(77)^ = 

(A;2-2)2A 

so that A = —1 and B = —A. Again considering Jacobi elliptic functions 
for which a, b are such that ^ — and j — are positive and also 



8 



D = ^2g~^ ; two solutions to (5.7) in terms of Jacobi elliptic functions (and 
their equivalent theta function representations) are 



(5.8) asciv) = J-T — 



A(2-P) \V2^' 



ll{2 - P) VV2^ 



A(20|-0^) ^^i(Wv/2^l^, 



where r is such that fc^ = 6'|/0|. 

Note that by the comments following equation (2.11), it is possible to ex- 
press the solutions obtained in this section in terms of Weierstrass functions. 
E. Abdalla and L. Correa-Borbonet [Ij have also considered the Einstein 
equation (i) with pif) = D/a(ty, and have found connections with Weier- 
strass functions in cosmic time (as opposed to the conformal time argument 
given here). In section 6 of this lecture we will find more general solutions 
to the conformal time equation (4.2) in terms of p, for arbitrary curvature k' 
and positive D-value. We will also consider the additional density functions 
p{t) = D/a{tf and p{t) = Di/a{tf + D2/a{tY for D, Di, D2 > 0. 
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6 FRLW and Weierstrass elliptic functions 



In general, if 51(0) = xq and g{u) satisfies 
(6.1) g\ur = F{g{u)) 

for F{x) — A/^x^ + A2,x^ + ^423;^ + Aix + Aq any quartic polynomial with no 
repeated roots, then the inverse function y{x) of g{u) is the elliptic integral 



(6.2) y{x) 



dt 



xo 



For initial condition g'{0) = 0, by (6.1) Xq is a root of the polynomial 
F[x). In this case the integral (6.2) can be rewritten as dzj P{z) where 
^ = l/{x — Xq) and P{z) is a cubic polynomial. To do this, first expand 
F{t) into its Taylor series about xq and then perform the change of variables 
z = l/(t — ,To). Fm^thermore one may obtain the form da/ \/Q{(t) for r = 
\{F' {xq)^+^F" {xq)) where Q{(j) = 4:a^ — g2<7—g3 is a Weierstrass cubic. This 
is done by setting z = (4(j-52/3)/B3 where B2 = F"{xo)/2 and B3 = F'{xo) 
are the quadratic and cubic coefficients oi P{z) respectively. Note that since 
F{x) has no repeated roots, xq is not a double root, F'{xo) ^ and the 

variable z is well-defined. Therefore we have obtained y ^ 24T-i"'(!:o) ^ ^ 
P~^{t). Writing this in terms of x, setting x — g{u) and solving for g{u), one 
obtains the solution to (6.1) 

(6.3) g{u) = xo + ^'^""'^ 



Mu; 92,93)-^ 



where 



1 . . 1 



2 



(6.4) g2 = AoA^ - -A.A^ + —A^ 

^3 = ^AoA2A^ + ^AiA2A3-^AlA^-^AoAl-^Al 

are referred to as the invariants of the quartic F{x). Since F{x) has no 
repeated roots, the discriminant A = g^f — 27g'| 7^ 0. Here if xq — 0, (6.3) 
becomes 

(6.5) g{u) = 



4^p{u; g2: gs) 
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For initial condition g'{0) ^ 0, by (6.1) xq is not a root of F{x) and 
a more general solution to (6.1) is due to Weierstrass. The proof (which we 
will not include here) was published by Biermann in 1865 (see OHHj). The 
solution is 



(6.6) 

v/FMp'(«) + |F'(xo) {p{u) - j-,F"{xo)) + ^F(xo)F'"(xo) 



g{u) =Xo + 



2 (piu) - j-,F"ixo)Y - ^F(xo)F(-)(xo) 



where p is formed with the invariants of the quartic seen in (6.4) such that 
A 7^ 0. Here if xq = 0, (6.6) becomes 

(6.7) g[u) = — 



2 {p{u) - ^A^Y - lAoA^ 

As a first example we consider the conformal time Einstein equation (4.2) 
for density p{t) = Di/a{tf + D2/a{tY with Di,D2 > 0. In this case (4.2) 
becomes 

(6.8) a'{rif = Aa(r/)^ - k'a{r]f + KdDia{ri) + K^Z^a 

and we take A^ = A, ^3 = 0, A2 = —k', Ai = KdDi and Aq = K^Da- 

The most general solution to (6.8) seen here is with initial conditions 

def. 

a'(0) 7^ so that a(0) =' ao is not a root of the polynomial 

(6.9) F{t) = At^ - k't^ + UdDit + KdD2. 

In this case the solution to (6.8) is given by (6.6) as 
(6.10) 

I , , ^^Wo)p\v) + \F'{a,) [piv) - l,F"{ao)) + l^F{a,)F"'{ao) 

2{p{v)-^,F"{ao)Y-i-,F{ao)F(^^Kao) 
with Weierstrass invariants 

(6.11) 92 = ATidD2 ' ^ ' 



12 



y-i 6 16 ^ 216 

restricted to be such that A = — 27g1 7^ so that F{t) does not have 
repeated roots. 



11 



Since D2, Hd > 0, by (6.9) zero is not a root of F{t) and therefore for 
initial conditions a'(0) 7^ and = the solution to (6.8) is given by (6.7), 

(6. 12) a(,) = v^PX'ri + r-f^ {pM + j|) 

2(p(,,) + i)'-lAS:,B2 

for invariants 5'2,5'3 as in (6.11) with A 7^ 0. One can compare this with the 
results in papers by Aurich, Steiner and Then, where curvature is taken to 
be k' = -1 [A, 5J. 

Finally, for a'(0) = and Oq 7^ a root of F{t) in (6.9), the solution to 
(6.8) is given by (6.3), 

(6-13) -i'l) = ao + 4p(,J'i'"°ik) 

again with invariants (6.11) such that A 7^ 0. 

For a more concrete example, consider the density function p{t) = D/a{t)^ 
for D > so that conformal time equation (4.2) becomes 

(6.14) a'{riY = Aa(r/)^ - k'a{riY + KdDa{r]). 

Here zero is a root of the polynomial F{t) with = A, ^3 = = 0, A2 = 
—k' and Ai = UdD. Therefore with initial conditions a'(0) = a(0) = 0, the 
solution to (6.14) is given by (6.5) as 

(6.15) a{r^) ^^'^ 



I2p{7]) + k' 
with invariants 

(6.16) 92=^^ and = -LIt^^d^ + ^ 
\ ) ^^12 16 216 

restricted to be such that A 7^ 0. 

As noted in the comments following equations (2.11), one can write this 
solution in terms of Jacobi elliptic functions (by using equations (2.11) if the 
roots of the reduced cubic 4t^ — g2t — gs are real). To demonstrate this, we 

choose D = 3^^-^/^ and fc' = 1 so that = and g2 = 1/12. For this 

positive curvature case, (6.15) becomes 

a/2 

(6.17) a{7])= ^ 



'3A(12p(r7) + l) 



12 



and the reduced cubic is At^ - {l/12)t = 4t(t- l/4v/3)(t+ l/4v/3). Applying 
(2.11) with eg = -1/4^3,62 = 0, ei = l/4v^, (6.17) can be equivalently 
written in terms of Jacobi eUiptic functions for rj > l/4-s/3 as 



(6.18) a{r]) 



-\/3 — 3 + 6ns^ 



2/A 



V3 + 6ds^ 



since fc^ = ~ I 7^ = ci — 63 = l/2-\/3. In terms of theta functions, 



(6.18) becomes 
(6.19) 

a{r]) 



{VS - 3) ei ej{rJ/^/2V3 91 r) + 6 9^ ^Kr^/v^ 91 r) 

^j2/lelej{r^/V2V3elT) 

VS 91 9l{r|/^/2V3 91 r) + 6 9^ 9l 9l{rj/^/27^ 9% r) 



2/A el el{r]/^2V^el,T) 



(v^ + 3) el eiir^ 1^271 el r) + 6 el ei{r^iVWl el r) 

where r is taken such that 1/2 = ^^2/^3- 

As a final example, we return to p(t) = D/a{ty for D > 0, which was 
considered in section 5. That is, we will obtain alternate solutions to equation 
(5.2). Since D > 0, zero is not a root of the polynomial F{t) with A4 — 
A,A^ = Ai = 0, A2 = —k' and Aq = KdD. Therefore for initial conditions 
a'(0) 7^ and a(0) = 0, (6.7) gives the solution 



(6.20) a{rj) = ^P'iv) 



13 



for invariants 



(6.21) g2 = KkdD + 



{k 



y\2 



12 



and 



y6 6 216 



restricted to be sucli that A 7^ 0. (6.20) is more general tlian tlie solutions 
to (5.2) in section 5, since here the curvature k' and the constant D are 
unspecified. 

To see this solution expressed in terms of Jacobi elliptic functions, take 
curvature k' = 1 and D = so that Qs = and g2 = 1/9. Then the 

reduced cubic is At^ - (l/9)t = 4(t - l/6){t + 1/6) so that 63 = -1/6,62 = 
0, ci = 1/6 and (6.20) becomes 

(6.22) a{r]) 



71 ''^"^ 



144((5(r,) + ^)2-l' 



Using (2.9) and (2.11) to write (6.22) in terms of Jacobi elliptic functions for 
77 > 1/6, 



(6.23) 



0(77) = 



2 - 3sn2 ( 4 



+ sn 



4 ( v 1 



n 1 



V •> V ^ 



where each of the positive and negative square roots solve (5.2) for k' — 
1,D^ l/36A/?d and where 7^ = ci-eg = 1/3 and modulus k'^ = = 1/2. 
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Writing (6.23) in terms of theta functions, 0(77) is 
(6.24) 



V6A ^2^4(77/ v^^i) {29j9l{ri/Vsei) - 9iei{rj / VseD) 



I29l9i9liv/y^0l) - 9l9l{v/V39l) 



2\/3A 9M{V/V39I) V 2919191(71/ V39I) + (t^/Vs^I) 
949392{v/V^9l)9,{r]/V39l) 



VQA^9l9f{ri/VS9l) + 39l9l9l{v/VS9l)9i{v/VS9l) + 29t9l{rj/VS9l) 

by the theta function representations for each of sn, ds and cs respectively. 
Here the dependence on r in the theta functions is suppressed and r is taken 
to satisfy 1/2 = 9l/9l 



A similar procedure can be done with (6.20) when k' — —1 and D 



1 



36AK(J 



SO that again gs — 0,g2 — 1/9 => 63 — —1/6, 62 = 0, ei = 1/6. In this case 
(6.20) is 



12 



(6.25) 



a{r)) 



144(p(r/)-^)2-l- 



In terms of Jacobi elliptic functions for 77 > 1/6, 
(6.26) 

'jL 



sn 



a{ri) 



6^-^ (*'7l)-3-^(*7s)+2 



1 


2ds^ 


' ri 


V2J 


+ 1 




2ds^ 


f rj 
i,73' 


72j 


- 1 



(*.75)+3-^(;S.75)+2 



sc 
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since again 7^ = ci — 63 = 1/3 and k'^ — — 1/2. Or, equivalently in 

terms of theta functions. 



0(77) 



(6.27) 

e2e'Mr]/V39l)9,ir]/V39l) 



e|e^iv/Vse|) 29',9l9l{r,/VS9l) + 9l9l{r,/VS9l) 



2VSA 9M{r,/V^9l) V 29l9l9l{rj / ^/^9l) - (^/v^^l) 

= ^0i9i{rj 1^/2,91) + Wl9l9l{r^/^/2>9l)9l{r^ 1^/2,91) + 2^1^1(7^/73^1) 

9Mril^0D 

Vel 9,92{r]/^/29l) {29l9l{r) 1^/291) + ^i^f (^/VS^D) 

where the suppressed r is taken to satisfy 1/2 = ^2/^3- 

In (6.20), one could also take D — and the invariants (6.21) become 

= 1 + and — ~k' -\- ^^j^. lik' — 1, the roots of the reduced cubic 
At^ - (13/12)^+ (35/216) are 63 = -7/12,62 = 1/6, ei = 5/12 and (6.20) 
becomes 



(6.28) a{ri) 



A (2(p(r7) + - \) 



which solves (5.2) in this case. Writing this solution in terms of Jacobi elliptic 
functions for 77 > 5/12, 

(6.29) 



^i'n) — ~ — 7=7 — 7 \ 7 

2VA (ns [ri, - sn [rj, j 
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2sd {rj, f ) 

^sc4(r?,^)+5sc2 (r^,f )+4 

since 7^ = ei — 63 = 1 and A;^ = ll-ll ~ ^/^- Writing these solutions 
equivalently in terms of theta functions, 



(6.30) 
a{rj) = 



2 Va 9Mv/9I) i9pliv/9i) - 9191{^/91)) 
2929,9l9,{v/9l)9s{v/9l) 

939,{r^/9l)y/9l9f{r^/9l) + 5^i^l^?(W^i)^i(W^i) + 4^l^t(r//^i) 



where the suppressed r satisfies 3/4 — ^2/^3- 

For D = and A-' = the invariants are 02 = 1 and 03 = and the 

roots of the cubic At^ — t are 63 = —1/2, 62 = 0, ei 1/2. In this case (6.20) 
becomes 

p'iv) 



(6.31) a{ri) 



a/a (2p(77)^ - I) 



which solves (5.2) in this case. Writing this solution in terms of Jacobi eUiptic 
functions for 77 > 1/2, 



(6.32) 



V2A (ns {v,j^)-sn (r;,^)) 
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since 7' 



ei - 63 = 1 and k'^ = 



1/2. In terms of theta functions, 



(6.33) 



V2A e^9,{^/9i) [eWMiQi) - eiej{r)/ei)) 

2929,9l9i{ri/9l)93{ri/9l) 
VA^i9l9t9l{r,/9l) - 9l9t{ri/9l) 



where r is taken to satisfy 1/2 = ^^|/^^|. 

For D = and k' = —1 the invariants are = 13/12 and = 35/216 

and the roots of the cubic At^ - (13/12)t - (35/216) are 63 = -5/12,62 = 
-1/6,61 = 7/12. In this case (6.20) is 



which solves (5.2) in this case. Writing this solution in terms of Jacobi elliptic 
functions for 77 > 7/12, 



(6.34) 



0(77) 




(6.35) 




2Va {ns {ri,l) + sn (r/, |)) 
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'2sd {ri, I) 

^3sc4(r7,i)+7sc2 (?7,|)+4 
2\/a(cs (77,|) (^'D) 
since 7^ = ei — 63 = 1 and k"^ = f^zf^ = 1/4. In terms of theta functions, 
(6.36) 

^^^^ ^ ^3gi(^/gi)Vgtg!(^/gi) - 59'20loi{v/oM{vm) + ^otoi(vm 

^ ^a^iC^/gj) V3g3^gt(r^/gi) + 7^1^1^^(77/^1)^1(77/^1) + ^BlOMW) 
where the suppressed r is such that 1/4 = ^|/^|. 
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7 Summary 



There are a number of ways to see that elhptic and theta functions solve 
the {d + 1)— dimensional Einstein gravitational field equations in a FRLW 
cosmology with a cosmological constant. Here we considered a scenario with 
no scalar field and with density functions p{t) = Di/a(t)^ + D2/a{tY for 
Di,D2 > 0, p(t) = D/a{t)^ and p(t) = D/a(t)'^ scaling in inverse proportion 
to the scale factor a{t). In these cases the first Einstein equation (i) takes the 
form a(t)^ = an expression containing negative powers of the cosmic scale 
factor a{t). At this point, one could have introduced the inverse function y{x) 
oia{t) to obtain an expression for y{x) as the integral of a power of x divided 
by the square root of a polynomial in x. That is, y{x) would be an elliptic 
integral that is not normal; other authors have taken this approach [H [16] . 

def. 

Here, we switched to conformal time by a change of variables a{ri) = a{f{r])). 
This produced an equation of the form a'{riY = an expression containing 
nonnegative powers of the conformal scale factor a{r]). 

After reviewing the definitions and properties of elliptic and theta func- 
tions in sections 2 and 3, we introduced the FRLW cosmological model in 
section 4. In section 5 for p{t) = D/a{tY, we obtained a differential equa- 
tion for a (77) containing only even powers of 0(77) and constructed solutions 
in terms of Jacobi elliptic functions, restricted to particular values of the 
constant D, parameterized by modulus < A; < 1. The equivalent theta 
function representations for these solutions were recorded, and we noted the 
special cases for which the elliptic solutions reduce to elementary functions 
and the corresponding solution in cosmic time was also computed. In sec- 
tion 6, we considered each of p{t) = Di/a{tY + D2/a(t)'^, p{t) = D/a{t)^ 
and p{t) = D/a{tY with various initial conditions and obtained solutions 
in terms of Weierstrass functions for general curvature k' and constants 
Di,D2,D > 0. By considering these solutions restricted to certain values 
of D > for p{t) = D/a(tY and p{t) = D/a{ty, we wrote 0(77) in these cases 
equivalently in terms of both Jacobi elliptic and theta functions. 

In current joint work with Floyd Williams [9], it is shown that elliptic 
functions also appear in the presence of a scalar field for both the 

FRLW and Bianchi I d— dimensional cosmological models with a nonzero 
cosmological constant and with a similar density function p{t) scaling in in- 
verse proportion to the cosmic scale factor a{t). There we note that the 
equations of each of these cosmological models can be rewritten in terms 
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of a generalized Ermakov-Milne-Pinney differential equation [T71 (TUl [H], a 
type which the square root of the second moment of the wave function of 
the Bose-Einstein condensate (BEC) also satisfies. From this we establish a 
direct mapping between cosmological and Bose-Einstein quantities. On the 
cosmological side of the FRLW-BEC correspondence, by imposing an equa- 
tion of state p^{t) = wpfj,{t) {w constant) on the density P0(t) and pressure 
p^{t) of the scalar field (f){t) one can obtain a differential equation involving 
elliptic functions on the BEC side of the mapping. 
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